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Illuminating the darkness 

More evidence from faint-star proper-motions for a cool and warm 
component to the local dark matter in the Galaxy * 
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Abstract. We present new evidence, based on faint HST proper-motions, for a bi-modal kinematic population of 
old white dwarfs, representative of the Thick-Disk and Halo of our Galaxy. This evidence supports the idea of 
a massive Halo comprised of faint and old white dwarfs, along with an extant population of Thick-Disk white 
dwarfs. We show how most of the required dark matter in the solar vicinity can be accounted for by the remnants 
from these two components together. 
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1. Introduction 



The claim by Oppenheimer et al. ( |2001 ) that they had 
found a significant population of Halo stars from a kine- 
matic survey toward the South Galactic Pole opened up 
an interesting discussion regarding the nature of dark mat- 
ter in the solar neighborhood. Their discovery seemed to 
corroborate earlier findings using very deep HST pho- 
tometry that also pointed out to the existence of an 
as yet unobserved component to the Halo, identified 
as very cool and ancient white dwarfs (Mendez and 
Minniti 2000). However, several authors have indicated 
that Oppenheimer's sample could be also interpreted as 
the tail of a 'warmer' white dwarf component (in the sense 
of having shorter cooling ages and therefore higher surface 
temperatures), better ascribed to the intermediate Thick- 
Disk population of the Galaxy (see e.g., Reid et al. 2001; 
Reyle et al. 2001). 

Obviously, disentangling the true nature of 
Oppenheimer's objects is important to understand 
their contribution to the dark matter problem in the 
Galaxy. In this paper we use HST data, deep photometry 
and proper-motions, to demonstrate that there is indeed 
a very important population of what appears to be 
Thick-Disk white dwarfs with a large mass density in the 
solar vicinity. These Thick-Disk objects could not have 
been identified before on the Hubble Deep Fields South 



and North due to their high Galactic latitude (\b\ > 50°), 
which prevented the appearance on these fields-of-view of 
an important population of these objects, characterized 
by a steep concentration toward the Galactic plane. 

A fundamental discovery presented in this paper is 
that, regardless of whether the objects found are Thick- 
Disk or Halo white dawrfs, they together account for most 
of the dark matter in the solar neighborhood. 

2. Description of the photometric and kinematic 
sample 

Our data is based on a proper-motion membership study 
of the nearest Globular Cluster NGC 6397 by King et 



al. ( [1998! KACP98 thereafter). They used HST proper- 
motions to segregate cluster members from field stars in 
order to define, characterize, and study the faint end of 
the main sequence on the cluster, in a seminal paper us- 
ing HST's WFPC accurate relative astrometry. Since their 
main science goal was to study the cluster, the field stars 
received comparatively little attention in that paper. In a 
previous paper they had also used the first-epoch expo- 
sures to study the cluster's white dwarf sequence (Cool et 
al. 1996). 



Send offprint requests to: Rene A. Mendez 

* Based on observations collected with the Hubble Space 
Telescope, which is operated by AURA, Inc., under contract 
with the National Science Foundation. 



More recently, Castellani et al. (2001) used King's data 
set to compare the star- and color-counts with Galactic 
models, constraining a number of large-scale Galactic pa- 
rameters, but no attention was given to the kinematic in- 
formation, as derived from the proper- motions. 

In this paper, we use the combined photometric and 
astromctric information to characterize, specifically, the 
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bluc-faint objects found in the field population of the 
KACP98 sample. The main motivation is the recent dis- 
covery of spectroscopically-confirmed old white dwarfs in- 
dicative of a rather massive Galactic component com- 
prised of these remnants, as previously suggested by the 
photometric study of the stellar content on the Hubble 
Deep Fields South and North by Mendez et al. Ql996p , 
and Mendez and Minniti ( J200C| , MMOO henceforth). 

Relative proper-motions for all field and cluster stars 
in a field close to the cluster NGC 6397 were derived by 
KACP98 using an earlier version of the methods described 
in Anderson and King (2000, and references therein). The 
final proper-motions for field stars are motions relative to 
a reference frame in which the cluster is at rest. Since the 
internal velocity dispersion for the cluster (< 5 km s^ 1 or 
0.4 mas y , where mas stands for milli-arcsec) is smaller 
than the proper-motions detected, at the distance of the 
cluster, the width of the cluster proper-motion provides a 
good indication of the uncertainties in the derived relative 
motions. 

The full sample was segregated by KACP98 into three 
sub-samples, based on a proper-motion and photometric 
selection criteria: NGC 6397 main-sequence stars (1,386 
stars), non-main sequence members (64 stars), and generic 
field stars (929 stars) on a 4.5 arcmin 2 field-of-view located 
at (l,b) = (338.1°, —12.0°), encompassing the range 15 < 
I < 25 and < V — I < 4.0. Here, we will concentrate our 
attention on the field stars of the KACP98 sample. 

The uniqueness of this data set is the availability of 
proper-motions for a a sample of very faint stars in a 
field with modest Galactic latitude, which then samples 
a significant path through the Thick-Disk of the Galaxy 
(Mendez & van Altena 1996). Another fortunate coinci- 
dence is that, unlike other studies being carried out with 
HSTcrn different Globular clusters, which focus on the in- 
ner core of the clusters, this field was acquired toward the 
outer envelope of NGC 6397, some 4.5 arcmin from the 
cluster center, far from the central 20 arcsec dominated 
by the cluster's power-law cusp, and therefore minimizing 
the 'contamination' from cluster stars. 



2.1. Photometric calibration 

KACP98 published photometry in the F555W and F814W 
ban ds of the 'flight-system', as defined by Holtzman et 
al. ( 1995 , Ho95). In order to compare our photometry 
with models, we need to convert this photometry to a 
standard system. However, because we only have magni- 
tudes in two bands we can not apply the iterative method 



used by Mendez et al (1996), or MMOO, to derive cali- 
brated photometry. Instead, in this paper, the calibration 
is done on a two-step process, by adjusting the zero-points 
of the flight-system (equation |l| below), then by apply- 
ing the flight-to-ground transformations (equation ^ which 
yields equation [| below), and then explicitly solving for a 
quadratic equation on the standard photometry. 



Specifically, the mags given by KACP98 are defined by 
(equation (7) of Ho95): 

WFPC2 = -2.5logiDN.s- 1 ) + Z FG + 2.5 log GR t (1) 

where WFPC2 is the flight magnitude, DN is the measured 
flux per unit time on the camera, Zfg are the zero-points 
as given in Table 6 of Ho95, and GRi is the detector's gain 
(slightly different for all three chips, hence the subscript, 
in our case this is ~7 e~ ADU -1 ). 

On the other hand, the flight-to-ground transformation 
is accomplished by equation (8) on Ho95, which we repeat 
here for completeness: 



SMAG = -2.h\og(DN s- 1 ) + Z FS + 2.h\ogGR l + 

Tlfs x SCOL + T 2 , FS x SCOL 2 (2) 

where SMAG and SCOL are the standard-system mag- 
nitudes and colors respectively, T\ and T 2 are linear 
and quadratic coefficients determined empirically which, 
along with the zero-points, are given in Table 7 of Ho95. 
Combining equations |l| and || above thus yields: 



SMAG = WFPC2 + {Z FS - Z FG ) + 2.5 log GR t + 

T ltFS x SCOL + T 2 , FS x SCOL 2 (3) 

In our case we have magnitudes in two passbands only, 
so that SMAG could be either V or I (the closest match 
to the F555W and F814W passbands, see Ho95). Using 
equation |^, and subtracting the terms for the I band from 
those of the V band, we end up with: 



V-I = SCOL = (^5551^-^8141^) + 

(Tl, F S,F5Q5W — Ti yF s, F 8l4w) X SCOL - 
(T2,FS,F555W ~ T2,FS,F8li\v) X SCOL 2 

which is a quadratic equation on SCOL. Solving for 
SCOL yields thus the standard color, which, combined 
with equation ^ yields also the individual apparent mag- 
nitudes. There is no iteration necessary, and the standard 
magnitudes are uniquely determined (this is not the case if 
we have more than two bands, in this case the best solution 
is obtained by the iterative process, which best satisfies (in 
a minimum-residual sense) equation ^| simultaneously, for 
all passbands). The corrections from the instrumental to 
the standard system are, in any case, quite small. For ex- 
ample for F14W= 15.87, and F555W-F814W= 1.03, one 
finds that 1= 15.83 and V-I= 1.043. Because these correc- 
tions arc quite small, errors in the transformation (which 
can be up to 0.1 mag, Ho95) are minimal (i.e., we could 
have as well used the instrumental magnitudes). The un- 
certainty in the final calibrated magnitudes is estimated 
to be less than 0.05 mag (Cool et al. 1996). 
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2.2. Proper-motions 

The proper-motions given by KACP98 are actually dis- 
placements on a 32-month baseline with respect to the 
mean displacement of the bona-fide cluster stars. These 
displacements were converted to actual relative proper- 
motions using the WFC plate-scale, and the baseline. This 
proper-motions were then converted to motions in RA and 
DEC by applying the position angle of the observations. 
As mentioned before, the motions are measured with re- 
spect to a system that has a zero mean motion for the clus- 
ter. In order to compute true absolute motions we need to 
correct for the proper- motion of the cluster, i.e.: 



Table 1. Properties of the bulk proper- motions in 
Galactic longitude and latitude for NGC 6397 main- 
sequence cluster members and field stars. 



^rel 



t^aba 



,l cl 

H'abs 



where yJ rel are the observed relative motions for field stars 
with respect to the cluster, fi^bs * s tne absolute motion of 
the cluster, and H abs is the absolute motion of the field 
stars. Fortunately, the motion for NGC 6397 is known, and 
it is given in the comprehensive compilation by Dinescu et 
al. ( |i999| ), which has n% aU cosS = 3.30±0.50 mas y" 1 , and 
fjbf = —15.20 ± 0.60 mas y _1 . We finally converted the 
motion to the Galactic system for an easier interpretation 
of the results. In this system, the mean absolute motion 
for the cluster would be [if = —11.46 mas y _1 , and 
$ = -10.52 mas y" 1 . 



3. Analysis of the sample 

The vector point diagram (VPD) for the whole sam- 
ple is shown in Fig. Fig. la includes all the data, 
whereas Fig. lb shows a detail of the inner core of the 
proper-motion distribution. The pluses indicate the clus- 
ter members, whereas the small black dots indicate the 
non-members. In addition, several other interesting ob- 
jects have been marked by the solid and open star sym- 
bols, as well as the filled square symbol. These objects will 
be discussed further below. This figure clearly shows that 
there is a good separation between cluster and field stars, 
although the boundary, specially in the interface between 
the two distributions, is somewhat blurry (fortunately this 
is of no consequence to our analysis) . 

Table ^ indicates the main properties of the proper- 
motion distributions shown in Fig. ^, split into Main- 
sequence cluster members, and non-members. The method 
to compute these parameters, accounting for outliers, has 
been described in Mendez et al. (2000). 

From Table | we can estimate the accuracy of the 
proper-motions directly from the dispersion in the cluster 
members (which should actually have a very small disper- 
sion if there were no observational errors). We obtain an 
overall uncertainty of about 2.0 mas y _1 in both coordi- 
nates. We also note that the mean motion for the cluster is 
what we would have expected from the absolute motion of 



Parameter MS cluster members Field stars 
Gal. Longitude mas y^ 1 mas y _1 


Mean 


-11.00 


-Z.80 


Mean error 


0.06 


0.14 


Dispersion 


2.06 


4.30 


Disp error 


0.06 


0.22 


Gal. Latitude 






Mean 


-10.78 


0.62 


Mean error 


0.05 


0.14 


Dispersion 


1.91 


4.04 


Disp error 


0.05 


0.20 



values are for the entire sample, there is a slight increase 
in proper-motion errors, as judged from the width of the 
cluster distribution, when going to fainter magnitudes, for 
example while the dispersion is 1.99 and 1.88 mas y _1 in /i/ 
and fib respectively for 16 < I < 22, these values increase 
to 2.22 and 2.66 mas y^ 1 respectively for 22 < I < 24. 

The VPD allows us to mark out some potentially 
interesting kinematic outliers for further consideration. 
We have drawn a more or less arbitrary section around 
the main bulk of the proper-motions for identification in 
the reduced proper-motion and color-magnitude diagrams 
(see below). There are, however, other interesting objects 
that fall within the outliers boundary, and these will be 
discussed later. 



the cluster (see Section 2.2), after correction from relative 
to absolute, and within the observational errors as given 
on the second row of the table above. While the above 



Recently, Oppenheimer et al. ( 2001 ) found a large pop- 
ulation of white dwarfs with kinematics representative of 
the Galactic Halo. These objects were pre-selected from 
a photographic proper-motion survey using UK Schmidt 
plates. They used the so-called 'reduced-proper-motion' 
(RPM) parameter H, used for the first time by Luyten 
many years ago. The interesting property of this param- 
eter is that while it is directly observable once we have 
proper-motions and photometry, it is mostly an indicator 
of luminosity for distinct stellar populations. This param- 
eter is defined as: 

H m =m + 51og/x + 5 = M + 5logV tan - blogK 

where m and /z are the apparent magnitude and (absolute) 
proper-motion, and M and Vt an are the absolute magni- 
tude (in the same band) and total tangential velocity for 
the star, while if is a constant. 

As it can be seen from its definition, different stellar 
populations (with characteristic luminosities and kinemat- 
ics) will occupy a different locus on, e.g., an H vs. color 
diagram, hence H can be used to distinguish distinct pop- 
ulations from its position on the RPM diagram. 

Fig. U shows the Hy vs. V-I diagram for our HST 
sample. The same objects found on the VPD have been 
marked on the RPM. Indeed, the RPM diagram allows us 
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[mas/yr] [mas/yr] 

Fig. 1. Vector point diagram for the whole sample (a), and for an enlargement (b), showing the most interesting 
objects found in the field population. The outliers boundary (solid-line polygon) has been drawn by eye. The objects 
outside the polygon are marked by an open circle in Figs. § and [| Pluses are for cluster stars, dots for field stars. 
Among these we have the extremely red objects (solid squares), and the blue faint population (filled-star symbol, six 
stars, and open star symbol, one star). Typical proper-motions uncertainties (~ 2 mas y _1 ) are marked by the large 
plus symbol in Fig. lb. 



to clearly distinguish some objects that did not stand-out 
in the VPD (e.g., the red stars or some of the faint blue 
objects), and vice- versa (the general kinematic outliers). 
The properties of the most interesting objects found (i.e., 
the extremely red stars, and the faint blue stars) are given 
in Table §. 

As it can be seen in Fig. ^, there is a group of six very 
red stars with 21 < / < 24 and colors V — I > 3.6. These 
objects are interesting because, with their faint magni- 
tudes and red-colors, they could be representatives of a 
very faint and nearby sample of low-luminosity late-type 
dwarfs or even brown-dwarfs. There is also a group of six 
faint blue stars (22 < / < 25, V - / < 1.0) with in- 
termediate kinematic, and one fainter blue object with 
even more extreme kinematics (labeled as 'Extreme' on 
Tabic These faint blue objects, as will be argued be- 
low, are part of two distinct populations, belonging to the 
Thick-Disk of the Galaxy and to the Halo respectively, and 
represent the late-stages in the evolution of the precursors 
of these stars when the Galaxy was formed. The proper- 
ties of all these objects (i.e., the extremely red stars, and 
the faint blue stars) are given in Table |[ 

An inspection of Table || indicates that the red stars 
are actually not very nearby: If they belong to the disk 
population (as predicted from all Galactic models, see e.g. 
Castellani et al. (2001)), then their kinematics (assuming 
a tangential velocity of 40 km s _1 , typical of this popula- 
tion) indicates distances larger than 600 pc, and absolute 
magnitudes brighter than Mj ~ 12.0, characteristic of M5 
to M6-type dwarfs (Gizis 1997; Leggett 1992, her Fig. 1). 
Their colors also indicate similar distances and luminosi- 



LO 
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V-I 



Fig. 2. Reduced proper-motion diagram for the field stars 
toward NGC 6397. Same symbols as for Fig. 0. Solid 
and dotted lines are the theoretical locus for 0.8 Mq 
white dwarfs of different cooling ages, from the models 
by Chabrier et al. ( 2000) ) for a tangential velocity of 60 
and 120 km s" 1 respectively. The six solid stars match 
well the theoretical expectation for a population typical 
of the Thick-Disk, whereas there is one object that seems 
to have Halo-like kinematics. Size of the errors is marked 
by the plus sign. 



ties. For example, adopting the color-absolute magnitude 



relationship given by Gizis (200C, his Table 4 and Fig. 11) 
or Leggett ( 1992j , her Fig. 15), we have for V — I ~ 4.0 a 
typical My ~ 16.5 or Mj ~ 12.5. We conclude that none 
of these red stars are very nearby nor particularly faint. 



Rene A. Mendez: Illuminating the darkness 



5 



Table 2. Properties of the red and faint blue stars in the 
field population. 
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Blue 


22.79 


0.60 


-18.24 


10.71 




23.28 


0.12 


-6.24 


-7.91 




24.15 


0.68 


9.71 


1.90 




24.34 


0.45 


8.12 


-2.57 




24.51 


0.98 


10.53 


5.26 




24.89 


0.75 


9.28 


-1.09 


Extreme 


24.03 


1.03 


-20.07 


-24.13 



On the other hand these objects can not be extragalactic, 
because they exhibit motions that are significantly dif- 
ferent from zero, within our estimated 2 mas y _1 error 
budget (see Table 0). 

The object with the most extreme-kinematic (the star 
far to the left on the VPD, see Fig. |) has a V-I= 2.85 
and 1= 19.56. The absolute magnitude for this color is 
predicted in the range 12.0 < My < 13.5, depending on 
metallicity (Gizis 1997). Thus, for a total proper- motion of 
74.43 mas y _1 one derives a tangential velocity on excess 
of 200 km s -1 , i.e. more than five times the expected tan- 
gential velocity of a typical Disk star. This object is thus 
quite interesting for follow-up spectroscopy as a potential 
Halo subdwarf (the same is true for the other kinematic 
outliers, 16 objects in total, located outside the polygon 
on Fig. [I]). Alternatively, if we insist on this object having 
a tangential velocity of 40 km s _1 , one derives a distance 
of about 120 pc, leading thus to an absolute magnitude of 
My ~ 17.0, i.e., more than three magnitudes fainter than 
normal Disk stars of the same color, an unlikely solution. 

The sample of faint blue objects represents an inter- 
esting population that deserves further attention, in view 
of the recent discovery by Oppenheimer et al. (2001) men- 
tioned before. Oppenheimer's study was limited to much 
brighter magnitudes (with R < 19.8) than our study 
(20 < I < 25, see Fig. ||), but they have spectral con- 
firmation. Our objects are much fainter, and it is unlikely 
that we will have any confirmatory spectroscopy on them 
until the large next-generation telescopes of 30 m or more 
become available. Table ^| shows that all these objects have 
motions that differ by more than 2.5 times from both the 
mean motion of the field or the cluster, even considering 
the large width of the field's proper-motion dispersion of 
more than 4 mas s _1 given in Table [j]. 

Even though we can not hope to have spectroscopy on 
these objects yet, we can use the RPM diagram, along with 
theoretical models to, at least, compare their position on 
this diagram with what we would expect for ancient white 
dwarfs. On Fig. || we have superimposed the latest theoret- 




V-I 

Fig. 3. Reduced proper-motion diagram for local samples 
of white dwarfs selected according to different criteria. 
Small dots are from the Legget, Ruiz and Bergeron ( 1998 ) 
sample with fx > 0.8 arc-sec y -1 . Plus and solid- square 
symbols are for the sample from Ruiz and Bergeron (2001 ) 



with fi > 0.2 arc-sec y and /i < 0.2 arc-sec y respec- 
tively. Open and solid stars is our sample (same as in 
Fig. |^), and triangles are the low-luminosty white dwarfs 
ESO 439-26 and CE 51 with extreme kinematics, as de- 
scribed in the text. Solid and dashed lines are the same as 
in Fig. |. 



ical cooling sequences for old white dwarfs from Chabrier 
et al. (2000). The continuous line is for a population with 
a mean tangential speed of 60 km s^ 1 , characteristic of the 
Thick-Disk (TD model), whereas the dashed line is for a 
speed of 120 km s^ 1 , typical of the Halo (Halo model). As 
can be seen the brightest blue stars fall very close to the 
expected location for old white dwarfs from the Thick- 
Disk. The faintest point however could, in principle, be 
fit by either model. However, we argue that this object is 
better fit by the Halo model rather than the TD model. 
Indeed, if this was a Thick-Disk star, then the correspond- 
ing (cooling) age (i.e., since it left the main-sequence) from 
the theoretical sequence (almost 14 Gyr), would imply 
that the Thick-Disk is quite older than what most mod- 
els predict. It will also imply an absolute magnitude for 
this object of My ~ 17.6, thus becoming the faintest de- 
generate objects ever found (as an example ESO 439-26, 
the faintest spectroscopically confirmed WD known has 
My ~ 17.4, Ruiz et al. 1995). Instead, if this object is 
placed on the Halo model, then its derived cooling age is 
11 Gyr, fully consistent with the age for the Galactic Halo, 
and with a reasonable luminosity of My ~ 15.8. We must 
emphasize here that the derived tangential velocities for 
these objects are used only as an indication of the popu- 
lation to which they belong, but these velocities are not 
used to determine distances directly (see Section 4). 
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3.1. Comparison to local and simulated samples 

Before we can proceed with any certainty on the proba- 
ble identification of the faint blue objects as belonging to 
the stellar Thick-Disk or Halo components of our Galaxy 
as described in the previous paragraph, we need to com- 
pare this sample of putative white dwarfs with the kine- 
matics of local samples of (Disk) white dwarfs. A sam- 
ple with well defined selection criteria, which has been 
used for a determination of the (Disk) white dwarf lumi- 
nosty function is provided by Legget, Ruiz, and Bergeron 
( |1998| , LRB98 thereafter). Fig. | shows the locus in the 
RPM diagram of the LRB stars in comparison with our 
sample. At first it would seem as if there is no difference 
between the two samples, thus implying that we are ac- 
tually looking at normal Disk stars. However, there is a 
big difference in the two samples: The LRB sample has 
a strong kinematic bias in favor of stars with quite large 
proper- motion (n > 0.8 arc-sec y -1 ), and thus the mean 
tangential velocity of the LRB89 sample is not representa- 
tive of the Disk population of white dwarfs as a whole, but 
rather it overestimates it. This then displaces the mean lo- 
cus of these objects on the RPM to larger values of Hy, 
making them appear similar to that of our sample (we 
must note that our sample has not been pre-selected ac- 
cording to a large proper- motion). The mean tangential 
velocity for the LRB8 sample is 113 km s" 1 with a dis- 
persion of 73 km s . It would of course be interesting 
to compare directly our sample with that of Disk white 
dwarfs that have not been selected according to kinemat- 
ics. Such a sample does exist (Holberg, Oswalt, and Sion 
2001, HOS01 thereafter), but unfortunately there are are 
no V-I colors available for these stars which would al- 
low us to make a direct comparison with our objects. 
From HOSOl's (volume- limited) complete sample of white 
dwarfs closer than 13 pc one finds that their mean tan- 
gential velocity is only 45 ± 5 (m.e.) km s _1 (this is ac- 
tually based on 43 out of their 49 stars closer than 13 pc 
due to the lack of proper-motions for a few of them) . We 
note here that volume-limited samples of main-sequence 
stars from Hipparcos give values close to 40 km s^ 1 
(Dehnen and Binney 1998), i.e., similar to those found 
for Holberg's sample (actually one would expect white 
dwarfs to have larger motions since they are older than 
the disk, and much of the velocity scattering is generated 
in the first Gyr (Wielen 1977). Indeed, numerical simula- 
tions using Reyle's model (private communication, model 
described in Reyle, Robin, and Creze 2001), for the same 
Galactic location as our field, and with the same magni- 
tude, color, and error constraints as those of our sample 
give a < Vtan >— 45.8 ± 1.9 (m.e.) km s _l , in excel- 
lent agreement with the results from HOS01 sample of 
nearby Disk white dwarfs. From this analysis we see that 
the mean tangential velocity for our possible Thick-Disk 
white dwarfs is more than 7a away from the mean of the 
local (simulated) sample. A further indication of the effect 
of kinematic selection is provided by the sample of Ruiz 



and Bergeron (2001, RB01 thereafter). They have pub- 



lished a sample of new nearby spectroscopically confirmed 
white dwarfs, selected from their own proper-motion sur- 
vey. For their sample with n > 0.2 arc-sec y _1 one finds 
< Vtan >= 60 ± 8 (m.e.) km s _1 , while for [i < 0.2 arc- 
sec y _1 one finds < V tan >= 43 ± 12 (m.e.) km s _1 (this 
last sample is incomplete in proper-motion, so this value 
is probably an upper limit). 

While comparing samples of putative disk white dwarfs 
to those presented in this paper, one also has to keep in 
mind that it is not unlikely that a few of the 'local' samples 
may also have representatives from the Thick-Disk and 
Halo (Fig |.). This is the case, for example of ES0439- 
26 (with Mi = 16.29 and V tan = 78 km s -1 , Ruiz et 
al. 1995) or CE51 (with Mj = 16.47 and V tan = 68 km s _1 , 
RB01). The mean tangential estimates derived from these 
local samples will be biased toward higher values due to 
these kinematic interlopers. Unfortunately, spectroscopy 
does not help here because all white dwarfs, irrespective 
of the metallicity of the parent population, look alike. 

We can finally compare the locus of our sample on 
the RPM diagram to that predicted by Galactic mod- 
els. Here we use the predictions from the model by Reyle 
et al. (2001, and private communication) based on the 
Besangon model of stellar populations in the Milky Way, 
which has been (succesfully) used to identify the kinemat- 
ics of Oppcnheimer's sample with that of Thick-Disk stars 
(see also Reid et al. 2001). Model simulations where run 
for the same Galactic location as that of our sample, with 
the same magnitude, color and proper-motion limits, as 
well as the same observational errors in all these quanti- 
ties. The field-of-view of the simulation was larger than 
that of our sample in order to have enough (simulated) 
stars in the comparison. Fig. |^ shows the locus for our 
stars, and those predicted by the model. As it is obvious 
from this figure, the location of our stars is totally in agree- 
ment with the predictions from the model for Thick-Disk 
stars. Disk stars on the other hand exhibit values with 
smaller tangential velocity (smaller Hv), as it has been 
discussed before. Of course all this evidence is, in a way, 
indirect (even if, as dicussed before, we had spectroscopy), 
and one could still insist that these stars are actually mem- 
bers of the Disk population. But, in this case, how can one 
explain that they obey such a distint kinematic behavior 
in the RPM diagram in comparison with the normal Disk 
stars? This obviously points to these objects as beeing part 
of a separate higher- velocity component. We thus conclude 
that, based on the evidence we have, these few blue faint 
stars are very likely true representatives of the Thick-Disk 
of our Galaxy. The implications for this will be discussed 
in the following Section. 

A final point is important when comparing the results 
from this paper with those found by MM00 for the HDF- 
S&N. Fig. H shows the calibrated color-magnitude diagram 
for our sample. As it can be seen, the faint blue objects 
are much less well defined, becoming mixed with the rest 
of the distribution, and making it difficult to distinguish 
between different kinematic components, a point to which 
we will come back in the next Section. The outliers on 
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Fig. 4. Reduced proper-motion diagram for simulated 
samples of white dwarfs selected according to the very 
same criteria as those of our sample. The simulatiosn are 
based on the Besancon Galactic model of stellar popu- 
lation synthesis (Celine Reyl/'e, private communication). 
Small dots are (simulated) Disk white dwarfs, plus sym- 
bols are (simulated) Thick-Disk stars, and star symbols 
are the objects found in this study. As it can been seen, 
our objects trace well the predicted locus for Thick-Disk 
stars as inferred from the model. 
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Fig. 5. Color-magnitude diagram for the whole sample 
of field stars, adopting the same symbols as for Figs. 1 
and 2. The separation between distinct kinematic groups 
becomes fuzzy. 



the VPD diagram (open circles) also become well mixed 
with the rest of the distribution, and only the red objects 
stand-out by virtue of their color-selection criteria. 

4. Discussion 

It is interesting to 'weight' (literally) the consequences of 
the presence of two kinematic components on the faint 
blue population described above. First of all, is the dis- 
covery of one possible truly Halo white dwarf (the open 
star on Figs. 1 and 2) consistent with previous findings? 



As it was shown by MMOO, the total mass contributed by 
a Halo-like population is given by: 



M{r max ) = < M WD > xN WD 
8.46 x l(T 8 ft x p Q 



r 2 ,!r ,41 



where pq is the local mass density of this component, 
f2 is the field-of-view in arcmin 2 , R is the galactocen- 
tric distance (in pc), Rq is the Solar Galactocentric dis- 
tance (assumed to be 8,500 pc), r is the heliocentric 
distance (pc), and M(r max ) is the mass (in Mq) con- 
tributed by this component out to a Heliocentric dis- 
tance r max (also in pc), while Nwd is the number of 
observed white dwarfs of mean mass < Mwd >■ If we 
compute the ratio of M(d max ) between, e.g., the Hubble 
Deep fields South and this field close to NGC 6397, then 
we find that, to the same depth in apparent magnitude, 
the expected number of Halo white dwarfs would scale as 
N WD (N 6397) ~ 1.11 X N W d{HDF - S), taken into ac- 
count the different fields-of-view of the two samples. Now, 
down to / ~ 25 we see from Fig. 3 on MMOO that we 
have two possible Halo white dwarfs on HDF-S, whereas 
we have found one based on its kinematic properties on 
the field of NGC 6397. No doubt, it is small numbers, but 
at least the two studies (MMOO based only on differential 
star counts, the present one on kinematics) are consistent 
with each other. Note also that the HDF images go about 
5 magnitudes fainter than the ones on the present study. 
We conclude thus that there is consistency between these 
results. 

In our study we seem to find a population of stars 
somewhat brighter than those detected by MMOO, with 
kinematics indicative of the Thick-Disk. However, there is 
no indication for these objects on the HDF-S&N sample 
by MMOO, which seems to be inconsistent with the find- 
ings of the present paper. Here we should point out that 
there was no kinematic information for the sample used 
by MMOO and, thus, it was not possible to kinematically 
distinguish Thick-Disk from Halo-like objects. Therefore, 
it is in principle possible that some of the blue objects 
on MMOO could actually be Thick-Disk stars. However, 
it turns out that the expected density of these latter ob- 
jects is much smaller in the HDF-S&N because of their 
spatial location in the Galaxy. If we adopt a density law 
characteristic of the Thick-Disk of the form: 



p( r ) = PQ exp" 



(5) 



where p(r) is the stellar densities at Galactic position r 
whereas pq is the density in the solar neighborhood, h z 
and h,R are the scale height and length of this compo- 
nent respectively, z is the distance from the Galactic plane, 
zq ~ 7(pc) is the Sun displacement with respect to the 
Galactic plane, and R p i is the distance from the Galactic 
center, as measured on the plane of the Galaxy, we find, 
using the analogous of equation |] for this density law, that 
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the ratio between the expected number of objects on the 
HDF-S&N and the field under consideration is only 0.29 
and 0.16 for HDF-S and HDF-N respectively. So, if we 
have 6 Thick-Disk stars on the N6397 field, we should only 
expect to see 1 or 2 such objects on HDF-S and 1 or none 
on HDF-N, for a reasonable choice of h z (760 pc) and hn 
(2500 pc) in the above equation. It is interesting to notice 
that, above V ~ 25 (i.e. the same magnitude range as our 
stars, see Fig. ||), there is indeed only one blue object on 
HDF-N and two blue objects on the HDF-S (see Fig. 3 on 
MM00). These objects have magnitudes and colors similar 
to those of our kinematically selected Thick-Disk sample 
in NGC 6397, and in agreement with the above numbers. 
However, lacking kinematics for the HDF sample, we can 
not advance much more in this analysis (see below). 

4.1. Density estimates 

We can use the six stars identified as members of the 
Thick-Disk to compute the stellar and mass density asso- 
ciated to this component. For this purpose we use the gen- 
eralized 1/Vmax method proposed by Tinney et al. ( 1993[ ). 
This method extends the classical 1/Knax to the case of 
non- (spatially) uniform samples, such as those found here. 
The method defines a (generalized) density-weighted vol- 
ume Igen as: 



v„ 



p(r) 

Vol PQ 



dVol 



(G) 



where p{r) and pq are the stellar densities at position r 
and in the solar neighborhood respectively, and the inte- 
gral is carried out over the volume available to the sample, 
given the observational constraints. 

Then, the total stellar density is given, just as in the 
classical method, by: 



PQ 



E 



Si 



gcn i 



where the factors fi account for the incompleteness of the 
sample (/, > 1), and the sum is carried over the individual 
data points. Meanwhile, the standard deviation on the 
derived overall density is given by: 



Similarly, if the sample is only complete to an upper 
proper-motion limit p u and a bright apparent magnitude 
m,b, the minimum distance for inclusion would be: 

r min = maxlp-^fi/n^p-Ho - 2 ^*-^] 

For our sample, the value for pi is determined by the ob- 
servational errors. As mentioned in Sect. 2.2, our obser- 
vational errors are on the order of 2 mas y~ x . Therefore, 
one could consider to have a firm detection for, say, three 
times the error, or 6 mas y _1 , which we will adopt as 
our value for /i/. Objects with motions larger than that 
would qualify within our selection criteria. Similarly, the 
fastest-moving object in the sample (see Fig[j].) has a total 
proper- motion of nearly 75 mas y _1 , value thus adopted 
for p u . It is very likely that the detection limit is larger 
than this, and it is determined by the epoch difference 
used to obtain the proper motions and the confusion limit. 
However, we have made tests and found that, even if one 
where to adopt r m i„ = 0, our results will change by less 
than 0.1%. The apparent magnitude limits were chosen 
as TOfc = 15 (see Fig. ||), and nif — 25 (see Table 1 on 
KACP98). 

For each of the six putative Thick-Disk stars identified 
in Table |^, we computed their absolute magnitudes in the 
I pass-band from their observed V-I color, by interpola- 
tion on Table 3 of Chabrier et al. ( pOOOp . The I and Mi 
were then used to obtain photometric parallaxes. We must 
emphasize here that the kinematics of the sample is not 
used at all to determine photometric distances or absolute 
magnitudes, rather the kinematics is used as an indication 
of belonging to a given population, which then allows us 
to adopt the correspondig density laws (e.g., equation M 
for calculating the effective sampling volume (equation |6j) . 
Therefore our results are not dependant on the exact value 
adopted for Vtan on either the Thick-Disk or the Halo. 

With the above values, the limiting distances r m j n and 
''max were computed, and the integral on equation ^| was 
evaluated by direct integration, using the adopted Thick- 
Disk density law, equation ^. In this case, equation ^ can 
be written as: 



Si 



The volume over which the integral in equation || is to 
be extended is determined by the minimum and maxi- 
mum distances (r m j n and r max respectively) over which 
any star can contribute to the sample (see, e.g., Mendez 
and Ruiz 2001). 

If we have a sample with a lower proper-motion limit 
pi and a faint apparent magnitude m/, r max is given by: 

r ma x - minlp-^p/p^p-HO - 2 ^-^) 

where p is the parallax, /i is the proper-motion and m the 
apparent magnitude. 



K™ = 8.46 x 10~ s O 



exp W" exp p W" r dr 



— l^nax x 5R 

Where V m & x is the classical maximum allowed volume 
for a magnitude and proper-motion selected sample (see, 
e.g., Mendez and Ruiz 2001), and 5ft is a (dimensionless) 
correction factor which accounts for the non-homogeneous 
distribution of stars in space (we shall see that, in any case, 
this correction factor is relatively small, close to one, for 
our sample). Both terms are given, respectively, by: 



V, 



max 



3.46 x 10 -8 x n x 



max 



(8) 



(7) 
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Table 3. Photometric distances and derived parameters 
for the faint blue stars on Table |[ 



Table 4. Adopted incompleteness fraction and derived 
values for 5ft. 



Objects 



Age 
Gyr 



Mi 
mag 



Dist 
pc 



pc 



^max 

pc 



V 

kms~ 



W Objedfe,„ 
kms~'Bhikms~ 1 



f 



Blue 



3.9 
1.0 
4.9 
2.6 
10.2 
5.8 



13.62 
12.44 
13.85 
13.22 
14.67 
14.06 



682 
1472 
1148 
1675 
929 
1466 



192 
198 
151 
190 
146 
183 



1888 
2474 
1698 
2270 
1164 
1542 



-59 
-43 
53 
64 
46 
64 



35 
-55 
10 
-20 
23 



69 
70 
54 
67 
51 



_1.085 
1.103 
1.166 
1.212 
1.290 
2.167 



0.9813 
0.9691 
0.9854 
0.9741 
0.9947 
0.9868 



1.3307 
1.4440 
1.2957 
1.4048 
1.1989 
1.2654 



ExtreiSfe 1.153 1.3206 1.3212 



Extreme 10.6 14.80 702 294 1097 



-67 



-80 



"TUT 



and 



5ft = 



£rr exp 



exp 



-{Rpl-RQl/hR r 2 



(9) 



The density law in equation [7] or equation || has two 
parameters, namely h z and Hr, which introduce some am- 
biguity in the computation of the integral. We have thus 
computed the integral (and the corresponding correction 
factor R) for several combinations of parameters in the 
range 760 < h z < 1300 pc and 2500 < h R < 3500 pc, 
and have chosen the extreme values as representative of 
the 'modeling errors' involved in the computation. Finally, 
the factors fi have been computed as the reciprocal of the 
completeness fraction as inferred from the artificial-star 
experiments reported by KACP98 (their Table 1). Table | 
summarizes all the adopted and derived values for the six 
Thick-Disk and one Halo stars from Table |^. The cooling 
age on the second column of Table || is the age (since be- 
coming a white dwarf) inferred from the same table on 
Chabrier et al. ( 2000) ) used for determining the absolute 
magnitudes. 

Because of the Galactic location of this field, it is 
easy to show that, with small errors, the (/i;,/^) proper- 
motions are projected almost directly into the V (along 
Galactic rotation) and W (towards the North Galactic 
Pole) velocity components respectively. These velocities, 
computed using the distance given in the fourth column 
of Table ||, are indicated, for reference, in the seventh 
and eigth columns of the same Table. Finally, in the last 
column we give the approximate tangential velocity (ne- 
glecting the small contribution from the U component of 
the velocity towards the Galactic center). As it can be 
seen from the Table, the mean tangential velocity is close 
to 65 kms^ 1 for the Thick-Disk stars, and more than 
100 kms^ 1 for the one possible Halo white dwarf, in agree- 
ment with the value found from Fig. [2| and larger than 
the value expected for normal nearby Disk white dwarfs, 
as already discussed in Section 3.1. 

Equation ^ has thus been integrated using the param- 
eters specified in Table |j| The adopted incompleteness 
factors and the extreme values derived for the integral 
in 5ft are given in Table |]. 5ft m i„ was derived adopting 
h z = 760 pc and ha = 3500 pc, whereas $imax was derived 
using h z = 1300 pc and hR = 2500 pc. In the case of the 



Halo star, instead of using a double-exponential density 
law, we have adopted an R~ 3 law similar to that used 
in MM00, which leads to equation [| above (to be more 
precise V gcn is given by equation |] divided by the local 
stellar density /?q, while the definition of 5ft is similar to 
that adopted for a double exponential. In this case how- 
ever, the free parameter in the density law is the axial 
ratio of the (flattened) Halo, which was varied in the in- 
terval 0.8 (oblate spheroid, 5ft m i„) to 1.0 (round spheroid, 

5ftmaa; ) ) • 

If (pQ min , cr pQ ) and {pQ max , o> Q ) are the den- 
sities and its error derived from R m in and R ma x respec- 
tively, then the final value for the total density and its 
error was calculated from: 



PQ 



'PS* 



,and, 



With these equations, and the values given in Table | we 
obtain, for the Thick-Disk white dwarf stars a value of: 



P(D\ 



(13.82 ±4.39) x 10~ 3 stars/pc 3 
(11.06 ±3.51) x 10~ 3 M & /pc 3 



For the one star that probably belongs to the Halo, the 
same computation above yields: 



PQwdh = (5.31 ± 3.50) x 10 -3 stors/pc 3 
= (4.25 ± 2.80) x 10~ 3 M Q /pc 3 

Obviously, we have much larger errors in this latter case 
due to the use of only one object. In both cases we have 
assumed a mean mass per white dwarf of 0.8 Mq. 

Because the correction factors 5ft given in Table |^ are 
all close to one, one might as well have computed the stel- 
lar densities using the classical 1/V max method, without 
correcting for the variation of stellar density as a function 
of distance away from the Sun. In this case, one does not 
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need to adopt any density law whatsoever, and the method 
is completely model independent. It is thus interesting to 
compare the previous values with the ones that one would 
have derived in this case. Dropping the correction factor 
3? in equation [j], one finds: 

Pq wdtd = (13.82 ± 4.39) x 1CT 3 stars/pc 3 
= (15.66 ± 3.51) x 10" 3 M Q /pc 3 

And for the Halo star: 

P Qwdh = (7.02 ± 6.54) x 1(T 3 stars/pc 3 
= (5.62 ± 5.23) x 10~ 3 M Q /pc 3 

As we see, the densities are slightly overestimated, 
but also the errors are larger. Still, we obtain basically 
the same results using either method. This is important, 
because it implies that, regardless of the details of the 
(model) density law adopted, the stellar and mass densi- 
ties are indeed quite large. The implications of these re- 
sults are discussed in the next Section. 

5. Conclusions 

The stellar and mass densities derived above imply that 
the Thick-Disk white dwarf stars, a population unac- 
counted for in the Hubble Deep Fields South and North 
due to their high Galactic latitude, contributes to about 
85% of the missing mass in the solar neighborhood, whose 
density in the solar neighborhood is estimated to be 
PQ DM -l-26xlCr 2 M /pc 3 .^ 

Additionally, the value derived for only one possibly 
Halo star, albeit having a large error, is similar to the value 
derived by MM00 from the larger sample of Halo white 
dwarfs found in the HDF-S&N. They found a density of 
Pq — 4.64 x 10 -3 Mq/pc 3 . Their value, however, has 
to be taken with caution because the lack of kinematic in- 
formation prevented a clear separation of the populations 
(compare Figs. || and || above). It has been mentioned al- 
ready that one might expect one or two Thick-Disk stars 
on the HDF-S and one or none in the HDF-N. If these 
objects were inadvertently identified by MM00 as Halo 
stars, then the estimated density for these objects on the 
Halo would have been too large. It is interesting to point 
out that on HDF-S there are two stars on the WD tracks 
that are quite a bit brighter than the rest of the blue-faint 
group. In the HDF-N there seems to be a grouping of faint 
WDs, and a brighter one, which could also be ascribed to 
the Thick-Disk (see Fig. 3 on MM00). In any case, proper- 
motions for both the HDF-S&N sample would help settle 
this uncertainty. 

If we combine the present results for Thick-Disk stars, 
with the earlier results from MM00 for Halo stars, the 
overall density on WD stars for both populations seems 
to be able to account for the whole of the (local) missing 
mass in the Galaxy. The existence of a massive component 



of Thick-Disk evolved white dwarfs has been proposed re- 
cently by Gates and Gyuk ( 2001 ) to help explain some 
of the (many) difficulties with a massive Halo of ancient 
stars. However, our results suggest that, in addition to 
this Thick-Disk component, there would also be a quite 
massive Halo of remnant stars. Even though their local 
densities are similar, the mass locked in each population is 
quite different, in account of their different stellar density 
functions (double exponential vs. power-law). The sim- 
ple integration of these density functions over the whole 
Galaxy lead to masses of: 

M\y £)-td = ^PQ WDTD h z h 3 R ~ 1O 9 M 

,and, 



M WD ^ H = Aitpq wdh R 3 q \n(R max /R r , 



9 x 10 10 M G 







for the Thick-Disk and Halo components respectively. 
Rmax ~ 20 kpc and R m m ~ 1 kpc in the last equation 
are the overall extension of the Halo and the Halo core- 
radius respectively. 

The total mass in these components is consistent with 
the mass required to produce the MACHO events, as indi- 
cated by Gates and Gyuk (2001), and is also not inconsis- 
tent with the total mass estimates for the Galaxy out to a 
distance of 20 kpc. We also note that with the small scale 
height of the Thick-Disk population (h z < 1.5 kpc), their 
contribution to the optical depth of microlensing events 
would be quite small, and thus the primary sources of the 
MACHO events would still be primordially the Halo white 
dwarfs. 

The possibility of a white dwarf dominated (dark) 
halo has been criticized from many fronts (see e.g. 
Richer 1999 and Gates and Gyuk 2001). Among the 
most important criticisms is that the precursor of 
these stars would have produced metals at a rate 
greater than observed (Gibson and Mould 1997) or that 
Galaxy halos at high redshift would be brighter than 
observed due to the the white dwarf precursors (Chariot 
and Silk 1995). While these are no doubt important 
issues, Chabrier (1999) and more recently Fontaine 
et al. (2001, see also Fontaine's C.S. Beals Lecture at 



http://www.astro.ubc.ca/WD-workshop/talks/indcx.html ) 
have shown that many of these points can either be over- 
come, or are doubtful criticisms of this scenario. One 
example is the chemical evolution problem: The yield of 
a zero-metallicity stars is largely uknown, nor has there 
been much modelling of the structure of these stars (e.g., 
a zero metallicity star never undergoes a helium flash, 
which must have an effect on the the yield of metals in 
the PNe phase). 

Even if the problems described above with this 
scenario persist, as pointed out by Lynden-Bell and 



Tout (2001) in the first 'Russell lecture' of the new mil- 



lennium, they could be avoided if the objects found are 
'pristine' white dwarfs, objects that never underwent nu- 
clear reactions nor followed the 'normal' path of stel- 
lar evolution. These objects (with masses in the range 
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0.2 < M/Mq < 1.1) could have formed gradually and 
very slowly by accretion onto planetary-sized precursor 
bodies in low-density regions such as cooling flows and 
galaxy halos. By radiating their energy before collapsing, 
these bodies would grow in mass resting on the zero- 
point energy of confined electrons following the uncer- 
tainty principle, and never getting a temperature high 
enough to start nuclear burning. The difficulty with this 
scenario is that, although theoretically possible, prelimi- 
nary results suggest that the required accretion rates are 
too slow to allow the formation of these objects in less 
than a Hubble time, which keeps the puzzle regarding the 
origin of these objects still open. 

As we have seen before, the assignement of the bulk 
of the faint blue stars found here to the Thick-Disk is 
actually irrelevant in the sense that one could as well 
have used the classical l/V maK method to derive space 
densities, without regard to the origin or association to a 
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